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Triangular  G2-Splines 


Hartmut  Prautzsch  and  Georg  Umlauf 


Abstract.  We  introduce  curvature  continuous  regular  free- form  surfaces 
with  triangular  control  nets.  These  surfaces  are  composed  of  quartic  box 
spline  surfaces,  and  are  piecewise  polynomial  multisided  patches  of  total 
degree  8 which  minimize  some  energy  integral.  The  Bezier  nets  can  be 
computed  efficiently  from  the  spline  control  net  by  some  fixed  masks,  i.e. 
matrix  multiplications. 


§1.  Introduction 

Most  methods  known  for  building  G^-free-form  surfaces  need  polynomials  of 
relatively  high  degree,  namely  0(k 2),  see  for  example  [2,3].  Only  recently  in 
1995  this  high  degree  was  beaten  by  two  methods  giving  Gfc-free-form  surfaces 
of  bidegree  2k  + 2 with  singular  [6]  and  regular  [4]  parametrizations,  respec- 
tively. These  low  degree  surfaces  can  be  represented  by  a control  net  [4]  or  a 
quasi  control  net  [6],  and  can  be  designed  so  as  to  allow  for  subdivision. 

In  this  paper  we  will  transfer  the  method  given  in  [4]  to  triangular  box 
splines.  Here  we  restrict  ourselves  to  G2-surfaces  which  are  the  most  important 
for  practical  applications  besides  G1-surfaces.  Further  details  and  the  general 
case  are  presented  in  [5,7]. 

This  paper  is  organized  as  follows.  In  Section  2 we  introduce  n-sided 
G2-patches.  These  patches  are  used  together  with  generalized  G2-box  spline 
surfaces  to  build  surfaces  of  arbitrary  topology.  How  the  free  parameters  in 
the  construction  can  be  used  to  generate  G2-splines  that  minimize  certain 
energy  functionals  and  how  these  G2-splines  can  be  generated  efficiently  will 
be  discussed  in  Section  3. 


§2.  p- Patches 

The  simplest  C2-box  splines  are  those  over  the  three-directional  grid  of  total 
polynomial  degree  four.  In  this  paper  we  consider  only  these  box  splines.  A 
quartic  box  spline  surfaces  has  a regular  triangular  control  net  and  each  of  its 
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Fig.  1.  Schematic  illustration  of  a quartic  box  spline  patch  (gray)  and  its  control 
net. 

polynomial  patches  is  determined  by  12  vertices  (called  control  points)  which 
are  arranged  as  in  Fig.  1. 

Furthermore,  we  can  identify  in  any  triangular  net  regular  subnets  of  the 
form  of  Fig.  1.  These  subnets  determine  patches  forming  a generalized  box 
spline  surface.  A generalized  box  spline  surface  has  holes  corresponding  to 
the  irregular  vertices  in  the  net.  An  example  is  shown  in  Fig.  2:  The  control 
net  (left)  contains  an  irregular  vertex  of  valence  8 and  the  generalized  box 
spline  surface  (right)  has  an  8-sided  hole. 


Fig.  2.  A triangular  net  with  a vertex  of  valence  8 (left)  and  the  corresponding 
generalized  quartic  box  spline  surface  with  an  8-sided  hole  (right). 


If  every  irregular  vertex  is  surrounded  by  at  least  three  rings  of  regular 
vertices,  every  irregular  vertex  corresponds  to  exactly  one  hole  in  the  general- 
ized quartic  box  spline  surface.  In  this  case  an  n-sided  hole  is  surrounded  by 
a complete  surface  ring  consisting  of  3n  patches. 

We  now  describe  how  to  fill  such  holes  with  regular  G2-surfaces: 

First:  for  any  n > 3,  n ^ 6,  we  define  a special  generalized  box  spline  surface 
that  lies  in  the  xy- plane  and  has  the  control  net  shown  in  Fig.  3 (left)  for 
n = 5.  Its  control  points  are  the  points 
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for  i = 1 ,n  and  j = 0, ...  ,3  and  k = 0, . . . , 3 — j,  where  Cj  = cos(2m/n ) 
and  Si  = sm(2m/n).  Thus  this  surface  consists  of  3n  patches,  say  x„+i,  . . . , 
X4„,  which  are  shown  schematically  in  Fig.  3 (right). 


Fig.  3.  The  control  net  of  xn+i,...,X4n  (left)  and  the  4 n planar  patches 
xi,...,X4„  (right). 

Second:  we  construct  n patches  xi, . . . , x„  filling  the  hole  left  by  the  patches 
xn+1, . . . , x4n , see  Fig.  3 (right).  Let 

x;(u,n,w)  = ]Tb  lijkBfjk{u,v,w) 

be  the  quartic  Bezier  representation  of  the  patch  xi,  where  u,  v,  w are  barycen- 
tric  coordinates  with  respect  to  some  reference  triangle,  i.e.  u>0,v>0,  w>0 
and  u + v + w = 1.  The  Bezier  points  of  x;  are  determined  such  that  Xf  has 
C2-contact  with  xn+i.  This  fixes,  say  b-jfc,  for  i = 0,1,2.  Further,  we  set 
bjjoo  = 0 and  b^*.  = b 2t2j,2k/^-  Fig-  4 shows  these  Bezier  points  for  n = 5. 
Note  that  the  scaling  differs  from  Fig.  3. 


Fig.  4.  The  Bezier  points  of  xj, . . . ,xn  for  n = 5. 

Lemma  1.  The  patches  xi, . . . ,X4„  are  regular  and  form  a surface  without 
self  intersections. 
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Ml  G-points 
H B-points 
I I R-points 


Fig.  5.  The  Bezier  points  of  a p-patch  for  n = 5. 

A proof  of  this  Lemma  can  be  found  in  [5] . 

Third:  for  any  polynomial  p : 1R2  —*  1R3  we  call  the  union  of  all  patches 

Pi(u,v,w)  = p(xi(u,v,w)),  i = l,...,4n, 

a p-patch.  In  the  sequel  we  only  consider  p-patches  of  degree  (4  or)  8 deter- 
mined by  a (linear  or)  quadratic  polynomial  p.  The  Bezier  points  of  such  a 
p-patch  are  illustrated  schematically  in  Fig.  5 for  n = 5. 

Since  for  i = n + 1, . . . , 2n  the  patch  x*  has  C2-contacts  with  x;_n,  x j+n 
and  Xj_j_2n  the  patch  p;  also  has  G2-contacts  with  pi_n,pi+n  and  Pi+2n-  Sim- 
ilarly, P2n+»  and  P3„+j+i  have  C2-contact  for  i = 1 ,...,n,  where  p4n+]  = 
P3n+i*  Moreover,  since  a p-patch  is  part  of  a polynomial  surface  each  pt, i = 
1, . . . ,n,  has  G2-contact  with  Pi_i,  where  po  :=  pn. 

The  Bezier  points  that  define  the  G2-conditions  between  the  patches 
Pi,. . . ,p„  are  marked  by  the  underlying  dark  area  in  Fig. 5.  We  call  them 
the  G-points  of  the  p-patch.  Leaving  these  points  fixed  and  changing  the 
other  Bezier  points  arbitrarily  such  that  all  C2-joints  between  adjacent  p* 
are  preserved,  we  obtain  a modified  p-patch.  In  general,  it  does  not  lie  on  a 
polynomial  surface,  but  we  will  still  call  such  a modified  p-patch  a p-patch. 

Theorem  2.  Any  n-sided  hole  of  a generalized  box  spline  surface  can  be  filled 
by  a p-patch  having  a C2-joint  with  the  generalized  box  spline  surface. 

Proof:  The  boundary  and  the  cross  boundary  derivatives  up  to  order  two  of 
an  n-sided  p-patch  are  determined  by  45n  Bezier  points.  We  call  these  the 
B-points  of  the  p-patch.  In  Fig.  5 they  are  marked  by  the  grey  area. 

The  B-points  can  be  changed  such  that  the  p-patch  fits  into  an  n-sided 
hole  of  a generalized  box  spline  surface  with  a G2-contact.  The  remaining 
points  without  the  G-points,  here  called  R-points,  can  then  be  adjusted  such 
that  any  patch  p„+i, . . . , P4„  of  the  p-patch  has  G2-contact  with  all  its  neigh- 
bours. Namely,  all  C2-conditions  involving  R-points  form  a linear  system 
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Fig.  6.  A possible  arrangement  of  the  A-,  B-,  D-  and  G-points  of  p,  U pj+n  U 
Pi+2n  G Pi+3nj  * = 1)  • • • > n- 

for  the  R-points.  The  matrix  of  this  system  is  square  if  we  add  enough  zero 
rows.  After  an  appropriate  permutation  of  its  columns  it  is  even  a block-cyclic 
matrix.  This  system  has  an  18n  parametric  solution.  Hence  there  are  18n 
R-points  that  can  be  chosen  arbitrarily.  We  call  them  A-points.  The  other  R- 
points  are  then  determined  by  the  A-,B-  and  G-points  via  the  C'2-constraints. 
We  call  these  the  D-points. 

Fig.  6 shows  a possible  choice  for  the  A-  and  D-points.  Note  that  this 
choice  is  not  unique.  □ 


§3.  Fair  p-Patches 

The  construction  of  a p-patch  that  fills  a hole  of  a generalized  box  spline 
surface  in  Theorem  2 is  such  that  different  coordinates  do  not  interfere  with 
each  other.  So,  without  loss  of  generality,  we  restrict  ourselves  to  scalar  valued 
p-patches  in  the  sequel.  Thus  a point  is  no  longer  a point  in  1R3,  but  in  H1. 

The  G-points  of  a p-patch  are  certain  Bezier  points  of  a reparametrized 
quadratic,  say 

2 2 — i 

i(x>y)  = J2'52qiix'yj- 

i=0  j= 0 

Hence  the  G-points  depend  linearly  on  the  six  coefficients  qij,  which  we  call 
the  Q-points. 

Further,  as  explained  in  the  proof  of  Theorem  2,  the  D-points  depend 
linearly  on  the  A-,  B-  and  G-points.  Thus,  if  we  consider  the  B-points  fixed, 
all  other  Bezier  points  of  the  p-patch  depend  linearly  on  the  six  values  q and 
the  18n  A-points. 

To  obtain  good  looking  surfaces  we  determine  these  6 + 18n  free  parame- 
ters such  that  the  p-patch  minimizes  a quadratic  fairness  functional.  We  tried 
several  functionals  including  the  thin  plate  energy  [1].  Judging  the  visual  ap- 
pearance of  the  surfaces  by  their  isophotes  we  got  the  best  results  with  the 
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Fig.  7.  An  initial  control  net  (left),  parameter  lines  of  the  resulting  G2-surface 
(middle),  top-view  of  the  surface  showing  isophotes  (right). 


functional 
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The  D-points  of  the  p-patch  depend  linearly  on  the  A-,  B-  and  Q-points. 
So  we  can  view  F as  a quadratic  functional  in  the  A-,  B-  and  Q-points. 

Since  F is  positive  definite,  it  is  minimal  for  fixed  B-points  if  its  deriva- 
tives with  respect  to  the  A-  and  Q-points  are  zero.  Differentiating  F = 0 with 
respect  to  the  A-  and  Q-points  leads  to  equations  that  are  linear  in  the  A-,  B- 
and  Q-points.  Solving  for  the  A-  and  Q-points  shows  that  the  Bezier  points 
of  the  p-patch  minimizing  F depend  linearly  on  the  B-points.  In  other  words, 
there  is  a matrix  Mn  depending  only  on  F and  n such  that  Mn b is  the  vector 
of  all  Bezier  points  if  b is  the  vector  of  all  B-points. 

Fig.  7 shows  an  example  for  the  G2-p-patch  construction.  The  initial 
triangular  control  net  has  an  irregular  vertex  of  valence  5.  The  isophotes 
confirm  that  the  resulting  surface  is  G2. 

Fig.  8 shows  a similar  example.  The  control  net  is  the  same  as  in  Fig.  7. 
However,  here  we  used  a p-patch  consisting  of  9 n,  n = 5,  rather  that  4n 
patches  to  fill  the  n-sided  hole  of  the  generalized  box  spline  surfaces. 


Fig.  8.  Parameter  lines  of  the  resulting  G2-surface  (left),  top- view  of  the  surface 
showing  isophotes  (right). 

A more  complex  example  is  shown  in  Fig.  9.  This  G2-surface  was  com- 
puted by  the  same  method  as  Fig.  7. 
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Fig.  9.  An  initial  control  net  (left),  the  generalized  quartic  box  spline  surface 
with  several  holes  (middle),  the  resulting  surface  where  every  hole  is 
filled  with  a p-patch  (right). 


Remark  3.  The  matrices  M3 , M4 , M5 , M7 , Ms  and  Mg  can  be  found  on  the 
website  http : / / i33www . ira . uka . de . 

Remark  4.  The  construction  above  can  be  generalized  for  generalized  box 
and  half  box  spline  surfaces  of  smoothness  order  2k  and  2k  — 1,  respectively, 
see  [5]. 
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